Sequence Alignment-Based Prediction of Myosin 7A: Structural Implications and Protein Interactions.

Myosin, a superfamily of motor proteins, obtain the energy they require for movement from ATP hydrolysis to perform various functions by binding to actin filaments. Extensive studies have clarified the diverse functions performed by the different isoforms of myosin. However, the unavailability of resolved structures has made it difficult to understand the way in which their mechanochemical cycle and structural diversity give rise to distinct functional properties. With this study, we seek to further our understanding of the structural organization of the myosin 7A motor domain by modeling the tertiary structure of myosin 7A based on its primary sequence. Multiple sequence alignment and a comparison of the models of different myosin isoforms and myosin 7A not only enabled us to identify highly conserved nucleotide binding sites but also to predict actin binding sites. In addition, the actomyosin-7A complex was predicted from the protein-protein interaction model, from which the core interface sites of actin and the myosin 7A motor domain were defined. Finally, sequence alignment and the comparison of models were used to suggest the possibility of a pliant region existing between the converter domain and lever arm of myosin 7A. The results of this study provide insights into the structure of myosin 7A that could serve as a framework for higher resolution studies in future.

Analysis of the Conformational Landscape of the N-Domains of the AAA ATPase p97: Disentangling the Continuous Conformational Variability in Partially Symmetrical Complexes.

Single-particle cryo-electron microscopy (cryo-EM) has been shown to be effective in defining the structure of macromolecules, including protein complexes. Complexes adopt different conformations and compositions to perform their biological functions. In cryo-EM, the protein complexes are observed in solution, enabling the recording of images of the protein in multiple conformations. Various methods exist for capturing the conformational variability through analysis of cryo-EM data. Here, we analyzed the conformational variability in the hexameric AAA + ATPase p97, a complex with a six-fold rotational symmetric core surrounded by six flexible N-domains. We compared the performance of discrete classification methods with our recently developed method, MDSPACE, which uses 3D-to-2D flexible fitting of an atomic structure to images based on molecular dynamics (MD) simulations. Our analysis detected a novel conformation adopted by approximately 2% of the particles in the dataset and determined that the N-domains of p97 sway by up to 60° around a central position. This study demonstrates the application of MDSPACE in analyzing the continuous conformational changes in partially symmetrical protein complexes, systems notoriously difficult to analyze due to the alignment errors caused by their partial symmetry.

Komagataella phaffii Erp41 is a protein disulfide isomerase with unprecedented disulfide bond catalyzing activity when coupled to glutathione.

In the methylotrophic yeast Komagataella phaffii, we identified an endoplasmic reticulum-resident protein disulfide isomerase (PDI) family member, Erp41, with a peculiar combination of active site motifs. Like fungal ERp38, it has two thioredoxin-like domains which contain active site motifs (a and a'), followed by an alpha-helical ERp29c C-terminal domain (c domain). However, while the a domain has a typical PDI-like active site motif (CGHC), the a' domain instead has CGYC, a glutaredoxin-like motif which confers to the protein an exceptional affinity for GSH/GSSG. This combination of active site motifs has so far been unreported in PDI-family members. Homology searches revealed ERp41 is present in the genome of some plants, fungal parasites, and a few nonconventional yeasts, among which are Komagataella spp. and Yarrowia lipolytica. These yeasts are both used for the production of secreted recombinant proteins. Here, we analyzed the activity of K. phaffii Erp41. We report that it is nonessential in K. phaffii, and that it can catalyze disulfide bond formation in partnership with the sulfhydryl oxidase Ero1 in vitro with higher turnover rates than the canonical PDI from K. phaffii, Pdi1, but slower activation times. We show how Erp41 has unusually fast glutathione-coupled oxidation activity and relate it to its unusual combination of active sites in its thioredoxin-like domains. We further describe how this determines its unusually efficient catalysis of dithiol oxidation in peptide and protein substrates.

An unusual dual sugar-binding lectin domain controls the substrate specificity of a mucin-type O-glycosyltransferase.

N-acetylgalactosaminyl-transferases (GalNAc-Ts) initiate mucin-type O-glycosylation, an abundant and complex posttranslational modification that regulates host-microbe interactions, tissue development, and metabolism. GalNAc-Ts contain a lectin domain consisting of three homologous repeats (α, ß, and γ), where α and ß can potentially interact with O-GalNAc on substrates to enhance activity toward a nearby acceptor Thr/Ser. The ubiquitous isoenzyme GalNAc-T1 modulates heart development, immunity, and SARS-CoV-2 infectivity, but its substrates are largely unknown. Here, we show that both α and ß in GalNAc-T1 uniquely orchestrate the O-glycosylation of various glycopeptide substrates. The α repeat directs O-glycosylation to acceptor sites carboxyl-terminal to an existing GalNAc, while the ß repeat directs O-glycosylation to amino-terminal sites. In addition, GalNAc-T1 incorporates α and ß into various substrate binding modes to cooperatively increase the specificity toward an acceptor site located between two existing O-glycans. Our studies highlight a unique mechanism by which dual lectin repeats expand substrate specificity and provide crucial information for identifying the biological substrates of GalNAc-T1.

Structural analysis of PTPN21 reveals a dominant-negative effect of the FERM domain on its phosphatase activity.

PTPN21 belongs to the four-point-one, ezrin, radixin, moesin (FERM) domain-containing protein tyrosine phosphatases (PTP) and plays important roles in cytoskeleton-associated cellular processes like cell adhesion, motility, and cargo transport. Because of the presence of a WPE loop instead of a WPD loop in the phosphatase domain, it is often considered to lack phosphatase activity. However, many of PTPN21's biological functions require its catalytic activity. To reconcile these findings, we have determined the structures of individual PTPN21 FERM, PTP domains, and a complex between FERM-PTP. Combined with biochemical analysis, we have found that PTPN21 PTP is weakly active and is autoinhibited by association with its FERM domain. Disruption of FERM-PTP interaction results in enhanced ERK activation. The oncogenic HPV18 E7 protein binds to PTP at the same location as PTPN21 FERM, indicating that it may act by displacing the FERM domain from PTP. Our results provide mechanistic insight into PTPN21 and benefit functional studies of PTPN21-mediated processes.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes.

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 µM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

Charged Amino Acids in the Transmembrane Helix Strongly Affect the Enzyme Activity of Aromatase.

Estrogens play critical roles in embryonic development, gonadal sex differentiation, behavior, and reproduction in vertebrates and in several human cancers. Estrogens are synthesized from testosterone and androstenedione by the endoplasmic reticulum membrane-bound P450 aromatase/cytochrome P450 oxidoreductase complex (CYP19/CPR). Here, we report the characterization of novel mammalian CYP19 isoforms encoded by CYP19 gene copies. These CYP19 isoforms are all defined by a combination of mutations in the N-terminal transmembrane helix (E42K, D43N) and in helix C of the catalytic domain (P146T, F147Y). The mutant CYP19 isoforms show increased androgen conversion due to the KN transmembrane helix. In addition, the TY substitutions in helix C result in a substrate preference for androstenedione. Our structural models suggest that CYP19 mutants may interact differently with the membrane (affecting substrate uptake) and with CPR (affecting electron transfer), providing structural clues for the catalytic differences.

Structure, dynamics, and stability of the smallest and most complex 71 protein knot.

Proteins can spontaneously tie a variety of intricate topological knots through twisting and threading of the polypeptide chains. Recently developed artificial intelligence algorithms have predicted several new classes of topological knotted proteins, but the predictions remain to be authenticated experimentally. Here, we showed by X-ray crystallography and solution-state NMR spectroscopy that Q9PR55, an 89-residue protein from Ureaplasma urealyticum, possesses a novel 71 knotted topology that is accurately predicted by AlphaFold 2, except for the flexible N terminus. Q9PR55 is monomeric in solution, making it the smallest and most complex knotted protein known to date. In addition to its exceptional chemical stability against urea-induced unfolding, Q9PR55 is remarkably robust to resist the mechanical unfolding-coupled proteolysis by a bacterial proteasome, ClpXP. Our results suggest that the mechanical resistance against pulling-induced unfolding is determined by the complexity of the knotted topology rather than the size of the molecule.

Functional and topological analysis of PSENEN, the fourth subunit of the γ-secretase complex.

The γ-secretase complexes are intramembrane cleaving proteases involved in the generation of the Aß peptides in Alzheimer's disease. The complex consists of four subunits, with Presenilin harboring the catalytic site. Here, we study the role of the smallest subunit, PSENEN or Presenilin enhancer 2, encoded by the gene Psenen, in vivo and in vitro. We find a profound Notch deficiency phenotype in Psenen-/- embryos confirming the essential role of PSENEN in the γ-secretase complex. We used Psenen-/- fibroblasts to explore the structure-function of PSENEN by the scanning cysteine accessibility method. Glycine 22 and proline 27, which border the membrane domains 1 and 2 of PSENEN, are involved in complex formation and stabilization of γ-secretase. The hairpin structured hydrophobic membrane domains 1 and 2 are exposed to a water-containing cavity in the complex, while transmembrane domain 3 is not water exposed. We finally demonstrate the essential role of PSENEN for the cleavage activity of the complex. PSENEN is more than a structural component of the γ-secretase complex and might contribute to the catalytic mechanism of the enzyme.

Dual-target inhibitors based on COX-2: a review from medicinal chemistry perspectives.

Inhibitors of COX-2 constitute a class of anti-inflammatory analgesics, showing potential against certain types of cancer. However, such inhibitors are associated with cardiovascular toxicity. Moreover, although single-target molecules possess specificity for particular targets, they often lead to poor safety, low efficacy and drug resistance due to compensatory mechanisms. A new generation of dual-target drugs that simultaneously inhibit COX-2 and another target is showing strong potential to treat cancer or reduce adverse cardiac effects. The present perspective focuses on the structure and functions of COX-2, and its role as a therapeutic target. It also explores the current state and future possibilities for dual-target strategies from a medicinal chemistry perspective.

The p97/VCP adaptor UBXD1 drives AAA+ remodeling and ring opening through multi-domain tethered interactions.

p97, also known as valosin-containing protein, is an essential cytosolic AAA+ (ATPases associated with diverse cellular activities) hexamer that unfolds substrate polypeptides to support protein homeostasis and macromolecular disassembly. Distinct sets of p97 adaptors guide cellular functions but their roles in direct control of the hexamer are unclear. The UBXD1 adaptor localizes with p97 in critical mitochondria and lysosome clearance pathways and contains multiple p97-interacting domains. Here we identify UBXD1 as a potent p97 ATPase inhibitor and report structures of intact human p97-UBXD1 complexes that reveal extensive UBXD1 contacts across p97 and an asymmetric remodeling of the hexamer. Conserved VIM, UBX and PUB domains tether adjacent protomers while a connecting strand forms an N-terminal domain lariat with a helix wedged at the interprotomer interface. An additional VIM-connecting helix binds along the second (D2) AAA+ domain. Together, these contacts split the hexamer into a ring-open conformation. Structures, mutagenesis and comparisons to other adaptors further reveal how adaptors containing conserved p97-remodeling motifs regulate p97 ATPase activity and structure.

Cryo-EM structure of the KLHL22 E3 ligase bound to an oligomeric metabolic enzyme.

CULLIN-RING ligases constitute the largest group of E3 ubiquitin ligases. While some CULLIN family members recruit adapters before engaging further with different substrate receptors, homo-dimeric BTB-Kelch family proteins combine adapter and substrate receptor into a single polypeptide for the CULLIN3 family. However, the entire structural assembly and molecular details have not been elucidated to date. Here, we present a cryo-EM structure of the CULLIN3RBX1 in complex with Kelch-like protein 22 (KLHL22) and a mitochondrial glutamate dehydrogenase complex I (GDH1) at 3.06 Å resolution. The structure adopts a W-shaped architecture formed by E3 ligase dimers. Three CULLIN3KLHL22-RBX1 dimers were found to be dynamically associated with a single GDH1 hexamer. CULLIN3KLHL22-RBX1 ligase mediated the polyubiquitination of GDH1 in vitro. Together, these results enabled the establishment of a structural model for understanding the complete assembly of BTB-Kelch proteins with CULLIN3 and how together they recognize oligomeric substrates and target them for ubiquitination.

Characterization of erythroferrone structural domains relevant to its iron-regulatory function.

Iron delivery to the plasma is closely coupled to erythropoiesis, the production of red blood cells, as this process consumes most of the circulating plasma iron. In response to hemorrhage and other erythropoietic stresses, increased erythropoietin stimulates the production of the hormone erythroferrone (ERFE) by erythrocyte precursors (erythroblasts) developing in erythropoietic tissues. ERFE acts on the liver to inhibit bone morphogenetic protein (BMP) signaling and thereby decrease hepcidin production. Decreased circulating hepcidin concentrations then allow the release of iron from stores and increase iron absorption from the diet. Guided by evolutionary analysis and Alphafold2 protein complex modeling, we used targeted ERFE mutations, deletions, and synthetic ERFE segments together with cell-based bioassays and surface plasmon resonance to probe the structural features required for bioactivity and BMP binding. We define the ERFE active domain and multiple structural features that act together to entrap BMP ligands. In particular, the hydrophobic helical segment 81 to 86 and specifically the highly conserved tryptophan W82 in the N-terminal region are essential for ERFE bioactivity and Alphafold2 modeling places W82 between two tryptophans in its ligands BMP2, BMP6, and the BMP2/6 heterodimer, an interaction similar to those that bind BMPs to their cognate receptors. Finally, we identify the cationic region 96-107 and the globular TNFα-like domain 186-354 as structural determinants of ERFE multimerization that increase the avidity of ERFE for BMP ligands. Collectively, our results provide further insight into the ERFE-mediated inhibition of BMP signaling in response to erythropoietic stress.

Deep generative models of LDLR protein structure to predict variant pathogenicity.

The complex structure and function of low density lipoprotein receptor (LDLR) makes classification of protein-coding missense variants challenging. Deep generative models, including Evolutionary model of Variant Effect (EVE), Evolutionary Scale Modeling (ESM), and AlphaFold 2 (AF2), have enabled significant progress in the prediction of protein structure and function. ESM and EVE directly estimate the likelihood of a variant sequence but are purely data-driven and challenging to interpret. AF2 predicts LDLR structures, but variant effects are explicitly modeled by estimating changes in stability. We tested the effectiveness of these models for predicting variant pathogenicity compared to established methods. AF2 produced two distinct conformations based on a novel hinge mechanism. Within ESM's hidden space, benign and pathogenic variants had different distributions. In EVE, these distributions were similar. EVE and ESM were comparable to Polyphen-2, SIFT, REVEL, and Primate AI for predicting binary classifications in ClinVar. However, they were more strongly correlated with experimental measures of LDL uptake. AF2 poorly performed in these tasks. Using the UK Biobank to compare association with clinical phenotypes, ESM and EVE were more strongly associated with serum LDL-C than Polyphen-2. ESM was able to identify variants with more extreme LDL-C levels than EVE and had a significantly stronger association with atherosclerotic cardiovascular disease. In conclusion, AF2 predicted LDLR structures do not accurately model variant pathogenicity. ESM and EVE are competitive with prior scoring methods for prediction based on binary classifications in ClinVar but are superior based on correlations with experimental assays and clinical phenotypes.

Gαs is dispensable for ß-arrestin coupling but dictates GRK selectivity and is predominant for gene expression regulation by ß2-adrenergic receptor.

ß-arrestins play a key role in G protein-coupled receptor (GPCR) internalization, trafficking, and signaling. Whether ß-arrestins act independently of G protein-mediated signaling has not been fully elucidated. Studies using genome-editing approaches revealed that whereas G proteins are essential for mitogen-activated protein kinase activation by GPCRs., ß-arrestins play a more prominent role in signal compartmentalization. However, in the absence of G proteins, GPCRs may not activate ß-arrestins, thereby limiting the ability to distinguish G protein from ß-arrestin-mediated signaling events. We used ß2-adrenergic receptor (ß2AR) and its ß2AR-C tail mutant expressed in human embryonic kidney 293 cells wildtype or CRISPR-Cas9 gene edited for Gαs, ß-arrestin1/2, or GPCR kinases 2/3/5/6 in combination with arrestin conformational sensors to elucidate the interplay between Gαs and ß-arrestins in controlling gene expression. We found that Gαs is not required for ß2AR and ß-arrestin conformational changes, ß-arrestin recruitment, and receptor internalization, but that Gαs dictates the GPCR kinase isoforms involved in ß-arrestin recruitment. By RNA-Seq analysis, we found that protein kinase A and mitogen-activated protein kinase gene signatures were activated by stimulation of ß2AR in wildtype and ß-arrestin1/2-KO cells but absent in Gαs-KO cells. These results were validated by re-expressing Gαs in the corresponding KO cells and silencing ß-arrestins in wildtype cells. These findings were extended to cellular systems expressing endogenous levels of ß2AR. Overall, our results support that Gs is essential for ß2AR-promoted protein kinase A and mitogen-activated protein kinase gene expression signatures, whereas ß-arrestins initiate signaling events modulating Gαs-driven nuclear transcriptional activity.

CaaX-motif-adjacent residues influence G protein gamma (Gγ) prenylation under suboptimal conditions.

Prenylation is an irreversible post-translational modification that supports membrane interactions of proteins involved in various cellular processes, including migration, proliferation, and survival. Dysregulation of prenylation contributes to multiple disorders, including cancers and vascular and neurodegenerative diseases. Prenyltransferases tether isoprenoid lipids to proteins via a thioether linkage during prenylation. Pharmacological inhibition of the lipid synthesis pathway by statins is a therapeutic approach to control hyperlipidemia. Building on our previous finding that statins inhibit membrane association of G protein γ (Gγ) in a subtype-dependent manner, we investigated the molecular reasoning for this differential inhibition. We examined the prenylation of carboxy-terminus (Ct) mutated Gγ in cells exposed to Fluvastatin and prenyl transferase inhibitors and monitored the subcellular localization of fluorescently tagged Gγ subunits and their mutants using live-cell confocal imaging. Reversible optogenetic unmasking-masking of Ct residues was used to probe their contribution to prenylation and membrane interactions of the prenylated proteins. Our findings suggest that specific Ct residues regulate membrane interactions of the Gγ polypeptide, statin sensitivity, and extent of prenylation. Our results also show a few hydrophobic and charged residues at the Ct are crucial determinants of a protein's prenylation ability, especially under suboptimal conditions. Given the cell and tissue-specific expression of different Gγ subtypes, our findings indicate a plausible mechanism allowing for statins to differentially perturb heterotrimeric G protein signaling in cells depending on their Gγ-subtype composition. Our results may also provide molecular reasoning for repurposing statins as Ras oncogene inhibitors and the failure of using prenyltransferase inhibitors in cancer treatment.

Fusion crystallization reveals the behavior of both the 1TEL crystallization chaperone and the TNK1 UBA domain.

Human thirty-eight-negative kinase-1 (TNK1) is implicated in cancer progression. The TNK1 ubiquitin-associated (UBA) domain binds polyubiquitin and plays a regulatory role in TNK1 activity and stability. No experimentally determined molecular structure of this unusual UBA domain is available. We fused the UBA domain to the 1TEL variant of the translocation ETS leukemia protein sterile alpha motif (TELSAM) crystallization chaperone and obtained crystals diffracting as far as 1.53 Å. GG and GSGG linkers allowed the UBA to reproducibly find a productive binding mode against its host 1TEL polymer and crystallize at protein concentrations as low as 0.2 mg/mL. Our studies support a mechanism of 1TEL fusion crystallization and show that 1TEL fusion crystals require fewer crystal contacts than traditional protein crystals. Modeling and experimental validation suggest the UBA domain may be selective for both the length and linkages of polyubiquitin chains.

Ubiquitin-specific protease 11 structure in complex with an engineered substrate mimetic reveals a molecular feature for deubiquitination selectivity.

Ubiquitin-specific proteases (USPs) are crucial for controlling cellular proteostasis and signaling pathways but how deubiquitination is selective remains poorly understood, in particular between paralogues. Here, we developed a fusion tag method by mining the Protein Data Bank and trapped USP11, a key regulator of DNA double-strand break repair, in complex with a novel engineered substrate mimetic. Together, this enabled structure determination of USP11 as a Michaelis-like complex that revealed key S1 and S1' binding site interactions with a substrate. Combined mutational, enzymatic, and binding experiments identified Met77 in linear diubiquitin as a significant residue that leads to substrate discrimination. We identified an aspartate "gatekeeper" residue in the S1' site of USP11 as a contributing feature for discriminating against linear diubiquitin. When mutated to a glycine, the corresponding residue in paralog USP15, USP11 acquired elevated activity toward linear diubiquitin in-gel shift assays, but not controls. The reverse mutation in USP15 confirmed that this position confers paralog-specific differences impacting diubiquitin cleavage rates. The results advance our understanding of the molecular basis for the higher selectivity of USP11 compared to USP15 and may aid targeted inhibitor development. Moreover, the reported carrier-based crystallization strategy may be applicable to other challenging targets.

Visualizing single-molecule conformational transition and binding dynamics of intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) play crucial roles in cellular processes and hold promise as drug targets. However, the dynamic nature of IDPs remains poorly understood. Here, we construct a single-molecule electrical nanocircuit based on silicon nanowire field-effect transistors (SiNW-FETs) and functionalize it with an individual disordered c-Myc bHLH-LZ domain to enable label-free, in situ, and long-term measurements at the single-molecule level. We use the device to study c-Myc interaction with Max and/or small molecule inhibitors. We observe the self-folding/unfolding process of c-Myc and reveal its interaction mechanism with Max and inhibitors through ultrasensitive real-time monitoring. We capture a relatively stable encounter intermediate ensemble of c-Myc during its transition from the unbound state to the fully folded state. The c-Myc/Max and c-Myc/inhibitor dissociation constants derived are consistent with other ensemble experiments. These proof-of-concept results provide an understanding of the IDP-binding/folding mechanism and represent a promising nanotechnology for IDP conformation/interaction studies and drug discovery.

Garcinolic Acid Distinguishes Between GACKIX Domains and Modulates Interaction Networks.

Natural products are often uniquely suited to modulate protein-protein interactions (PPIs) due to their architectural and functional group complexity relative to synthetic molecules. Here we demonstrate that the natural product garcinolic acid allosterically blocks the CBP/p300 KIX PPI network and displays excellent selectivity over related GACKIX motifs. It does so via a strong interaction (KD 1 µM) with a non-canonical binding site containing a structurally dynamic loop in CBP/p300 KIX. Garcinolic acid engages full-length CBP in the context of the proteome and in doing so effectively inhibits KIX-dependent transcription in a leukemia model. As the most potent small-molecule KIX inhibitor yet reported, garcinolic acid represents an important step forward in the therapeutic targeting of CBP/p300.